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ABSTRACT 

The internal properties of the neutron star crust can be probed by observing 
the epoch of thermal relaxation. After the supernova explosion, powerful neutrino 
emission quickly cools the stellar core, while the crust stays hot. The cooling wave then 
propagates through the crust, due to its finite thermal conductivity. When the cooling 
wave reaches the surface (age 10 — 100 yr), the effective temperature drops sharply from 
250 eV to 30 or 100 eV, depending on the cooling model. The crust relaxation time is 
sensitive to the (poorly known) microscopic properties of matter of subnuclear density, 
such as the heat capacity, thermal conductivity, and superfluidity of free neutrons. We 
calculate the cooling models with the new values of the electron thermal conductivity 
in the inner crust, based on a realistic treatment of the shapes of atomic nuclei. 
Superfluid effects may shorten the relaxation time by a factor of 4. The comparison of 
theoretical cooling curves with observations provides a potentially powerful method 
of studying the properties of the neutron superfluid and highly unusual atomic nuclei 
in the inner crust. 
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1 INTRODUCTION 

Neutron stars are natural astro-physical laboratories of su- 
perdense matter. In their cores, at densities above the nu- 
clear matter density po = 2.8 x 10 14 g cm -3 , the properties 
of matter such as equation of state and even composition 
are largely unknown (e.g., Lattimer & Prakash 2000, and 
references therein). In the absence of exact theory of super- 
dense matter, different theoretical models predict different 
equations of state and compositions (neutrons, protons and 
electrons; hyperons; pion or kaon condensates; deconfined 
quarks) . 

One of the potentially powerful methods to probe the 
internal structure of isolated neutron stars is modelling of 
their cooling (e.g., Pethick 1992; Page & Applegate 1992; 
Page 1998a,b). The theoretical cooling curves depend on the 
adopted physical models of the stellar interior, especially the 
neutrino emission and heat capacity, as well as the superflu- 
idity of neutrons and protons in the core. Confronting the- 
ory and observations allows one, for example, to constrain 
the range of critical temperatures of the superfluidity (e.g., 
Yakovlev et al. 1999). 

Observing thermal emission of the very young neutron 
stars, t £ 100 yr, opens a possibility of studying the proper- 
ties of the neutron star crust. Soon after a supernova explo- 
sion, the young star has large temperature gradients in the 
inner parts of the crust. While the powerful neutrino emis- 



sion quickly cools the core, the crust stays hot. The heat 
gradually flows inward on a conduction timescale, and the 
whole process can be thought of as a cooling wave propaga- 
tion from the center towards the surface. During this thermal 
relaxation the effective temperature stays almost constant at 
about 250 eV. When the cooling wave reaches the surface, 
the effective temperature drops sharply by as much as an 
order of magnitude in the fast cooling scenario, and by a 
factor of 2 — 3 in the slow cooling scenario. The duration of 
the relaxation epoch depends mainly on the heat capacity 
and thermal conductivity of the inner crust (Lattimer et al. 
1994). 

Although the equation of state of matter at subnuclear 
density is known sufficiently accurately (Negele & Vautherin 
1973, Pethick & Ravenhall 1995), the properties of atomic 
nuclei are not. The nuclei become unusually neutron-rich, 
with the smooth proton and neutron distributions. At the 
bottom of the crust, p <; 10 14 g cm~ 3 , the nuclei can be 
nonspherical and form clusters (Lorenz et al. 1993, Pethick 
& Ravenhall 1995) . The liquid of neutrons dripped from the 
nuclei may be superfluid, with the critical temperatures that 
are very model-dependent. The thermodynamic and trans- 
port properties of this matter are subject of large theoretical 
uncertainty. 

In this paper we make improved models of the young 
neutron stars. We obtain the new values of the electron 
thermal conductivity in the inner crust, based on a realistic 
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treatment of the shapes of atomic nuclei. Using a new nu- 
merical code, we calculate the cooling models and determine 
the duration of the thermal relaxation epoch. We extend 
the analysis of Lattimer et al. (1994) and derive the depen- 
dence of the relaxation time on the microscopic parameters 
of the crust. Our preliminary results have been summarized 
by Yakovlev et al. (2000). 



2 COOLING MODEL 

2.1 Equations of thermal evolution 

Neutron stars are born very hot in supernova explosions, 
with the internal temperature T ~ 10 11 K, but gradually 
cool down. About twenty seconds after the birth, they be- 
come fully transparent for the neutrinos generated in numer- 
ous reactions in stellar interiors. We consider the cooling in 
the following neutrino-transparent stage. The cooling is re- 
alized via two channels, by neutrino emission from the entire 
stellar body and by heat conduction from the internal lay- 
ers to the surface resulting in thermal emission of photons. 
For simplicity, we neglect the possible reheating mechanisms 
(frictional dissipation of the rotational energy, Ohmic decay 
of the internal magnetic field, or the dissipation associated 
with weak deviations from the chemical equilibrium; see, 
e.g., Page 1998a). 

The internal structure of neutron stars can be re- 
garded as temperature- independent (e.g., Shapiro & Teukol- 
sky 1983). The relativistic equations of thermal evolution 
include the energy and flux equations (Thorne 1977): 



Table 1. Neutron star models 
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where Q v is the neutrino emissivity [ erg cm s J , Cv is 
the specific heat capacity [erg cm" 3 K _1 ], k is the thermal 
conductivity, and L r is the "local luminosity" defined as 
the non-neutrino heat flux transported through a sphere of 
radius r. The gravitational mass m(r) and the metric func- 
tion <!>(r) are determined by the stellar model. After thermal 
relaxation, the redshifted temperature T(t) = T(r, i) e*' 1 "- 1 
becomes constant throughout the interior. 

It is conventional (Gudmundsson et al. 1983) to subdi- 
vide the calculation of heat transport in the neutron star 
interior (r < Rb) and in the outer heat-blanketing envelope 
(Rb < r < R), where R is the stellar radius, and the bound- 
ary radius Rb corresponds to the density pb = 10 10 g cm" 3 
(~ 100 meters under the surface). The thermal structure 
of the blanketing envelope is studied separately in the sta- 
tionary, plane-parallel approximation to relate the effective 
surface temperature T 3 to the temperature Tj, at the inner 
boundary of the envelope. We use the T s -Tb relation ob- 
tained by Potekhin et al. (1997) for the envelope composed 
mostly of iron. 

The effective temperature determines the photon lu- 
minosity: L 7 = L r (R,t) = 4-7vaR 2 Tf(t). Both L 7 and T s 
refer to the locally-flat reference frame on the surface. A 
distant observer would register the "apparent" luminosity 
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a threshold configuration for switching on direct Urea process 
b configuration with maximum allowable mass 
c central density in 10 14 g cm -3 
d AR C rust is defined as R — Rcorc 



= L~f (1 — r g /R) and the "apparent" effective tempera- 
ture T s °° = T s y/l-r g /R, where r g = 2GM/c 2 is the grav- 
itational radius. Typically, r g /R — 30 — 40%. 

We have developed a new evolutionary code, based on 
the Henyey-type scheme on a grid of spherical shells (Kip- 
penhahn et al. 1967) . The hydrostatic model of the neutron 
star with a given equation of state is calculated separately 
and is fixed throughout the calculation. In the initial con- 
figuration the star has a constant redshifted temperature 
throughout the interior, T = 10 10 K, and no heat flux, L r — 
0. Also, to improve numerical convergence the thermal con- 
ductivity in the core is boosted for the initial epoch t < 10 -2 
yr. Since the crust is thermally detached from the core at 
such small age, this correction has no effect on the cooling 
curves. Full details of the new code are available on the in- 



ternet (http:/ /www. ast.cam.ac.uk/~ognedin/ns/ns. html) 



2.2 Physics input 

We use the equation of state of Negele & Vautherin (1973) 
in the stellar crust with the smooth composition model of 
ground-state matter to describe the properties of atomic nu- 
clei (Kaminker et al. 1999). We assume that the nuclei are 
spherical throughout the entire crust. The core-crust inter- 
face is placed at p cc = 1.5 x 10 14 g cm" 3 . For simplicity, we 
consider the neutron star cores composed of neutrons (n), 
protons (p) and electrons (e) and use the moderately stiff 
phenomenological equation of state of Prakash et al. (1988) 
(in the simplified version proposed by Page & Applegate 
1992), in agreement with our previous work (Yakovlev et al. 
1999 and references therein). 

The parameters of the models are summarized in Table 
[l]. It shows the stellar masses, radii, central densities, crust 
masses, and crust thicknesses for a number of models. The 
maximum mass of the stable neutron stars with this equa- 
tion of state is 1.73 Mq. We define the crust thickness as 
AJ? crus t = R — -Rcore, while the pr oper geom etrical thickness 
(for AR cruBt < R) is AR crust / y/l — r g /R. As we increase 
the mass M, the radii and crust masses of the stable config- 
urations get smaller, i.e., the stars become more compact. 

Free neutrons in the inner crust and both neutrons and 
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Figure 1. Density dependence of the critical temperatures of 
superfluidity (SF) of free neutrons in the inner crust, and neu- 
trons and protons in the core for the strong (solid lines) and 
weak (dashed lines) superfluid models (see text for details). 
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protons in the core of a neutron star are likely to be super- 
fluid. We assume the singlet-state pairing of the protons in 
the core. Superfluidity of free neutrons in the crust and of 
neutrons in the outermost part of the core is known to be 
of the singlet-state type, but at higher densities it switches 
to the triplet-state type. Various microscopic theories pre- 
dict a large scatter of the critical temperatures of the neu- 
tron and proton superfluids, T cn and T cp , depending on the 
nucleon-nucleon potential model and the many-body theory 
employed (see Yakovlev et al. 1999 for references). 

As an example we will use two models, a weak and a 
strong superfluidity (Fig. |l|). The model of strong superflu- 
idity corresponds to the higher T c . It is based on the rather 
large energy gaps calculated by Elgar0y et al. (1996) for 
the singlet-state pairing (with the maximum gap of about 
2.5 MeV) and by Hoffberg et al. (1970) for the triplet-state 
pairing. The weak superfluid model makes use of the smaller 
gaps derived by Wambach et al. (1993) (with the maximum 
gap of about 1 MeV) for the singlet-state superfluid and by 
Amundsen & 0stgaard (1985) for the triplet-state neutron 
superfluid. For simplicity, we use the same function T c (n) 
to describe the singlet pairing of free neutrons in the crust 
(n = n n ) and of the protons in the core (n = n p ). 

The cooling of neutron stars is mainly determined by 
the neutrino emissivity, specific heat capacity, and thermal 
conductivity. We include all relevant sources of neutrino 
emission (Yakovlev et al. 1999, 2001): the direct and mod- 
ified Urea processes, nn, pp and np bremsstrahlung in the 
core; and plasmon decay, e~e + pair annihilation, electron- 
nucleus (eZ) and nn bremsstrahlung in the crust. The emis- 
sivity of the proton branch of the modified Urea process 
derived by Yakovlev & Levenfish (1995) has been corrected 
by Yakovlev et al. (2001), which has almost no effect on 
the cooling curves. We include the proper reduction of the 
neutrino reactions by superfluidity and also an additional 
neutrino emission due to Cooper pairing of superfluid nucle- 
ons. The effective nucleon masses in the core and the crust 
are set equal to 0.7 of their bare masses. 

A very important cooling effect is produced by the pow- 



Figure 2. Density dependence of the neutrino emissivity at 
T = 10 9 K. Solid lines: partial emissivities due to electron-nucleus 
bremsstrahlung (e-brems) and plasmon decay (plasma) in the 
crust and the total emissivity produced by direct and modified 
Urea processes and by nucleon-nucleon (NN) bremsstrahlung in 
the nonsuperfluid core. Dashed and dash-and-dot lines: the emis- 
sivity due to Cooper pairing of neutrons for the models of strong 
and weak superfluidity in the crust. Vertical dotted lines indicate 
the neutron drip density and the boundary of the core. 
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Figure 3. Density dependence of the specific heat capacity at 
T = 10 9 K. Solid lines: partial heat capacities of ions (i), elec- 
trons (e) and free neutrons (n) in nonsuperfluid crusts, as well 
as of neutrons, protons (p) and electrons in nonsuperfluid cores. 
Dashed lines: heat capacities of free neutrons in the crust modified 
by weak or strong superfluidity. 
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Figure 4. Density dependence of the electron thermal conduc- 
tivity at T = 10 7 , 10 s and 10 9 K in the core and the crust. 

erful direct Urea process (Lattimer et al. 1991). For our 
equation of state, this process is allowed at densities above 
the threshold, p cr it = 1.298 x 10 15 g cm -3 . If the central den- 
sity of the model exceeds the threshold, p c > p cr it , the stellar 
core has a central kernel where the direct Urea process leads 
to fast cooling. The masses and radii of these kernels, Mo 
and 7?d , are given in Table ^. The mass of the central kernel 
increases rapidly with M. In addition, we show the thresh- 
old configuration with p c — pcrit • It separates the low-mass 
models, where the direct Urea process is forbidden, from the 
high-mass models, where the direct Urea is allowed. 

For illustration, Fig. ^ shows the emissivity of various 
neutrino processes versus density at T = 10 9 K. The plas- 
mon decay is a powerful neutrino emission mechanism in 
a hot crust at not very high densities, but it fades away 
quickly when the temperature decreases below 10 9 K. The 
electron-nucleus bremsstrahlung is efficient throughout the 
entire crust. In superfluid crusts, the Cooper pairing neu- 
trino emission switches on at temperatures T = T cn , reaches 
maximum at T slightly below T cn , and fades away exponen- 
tially for T <C Ten- In the neutron star core, a large jump 
of the neutrino emissivity at p = p cr it is associated with the 
direct Urea process. 

The heat capacity is contributed by neutrons, protons 
and electrons in the core; and by electrons, free neutrons, 
and atomic nuclei (vibrations of ions in Coulomb lattice) 
in the crust. The superfluid effects on the heat capacity of 
nucleons in the core and of free neutrons in the crust are 
incorporated according to Levenfish & Yakovlev (1994). In 
the absence of superfluidity, the neutrons would have the 
dominant contribution in the core and the inner crust. The 
effects of neutron superfluidity are illustrated on Fig. ^| in 
the case of crust superfluidity. When T falls only slightly 
below T cn , the superfluidity increases the neutron heat ca- 
pacity due to the latent heat released at the phase transi- 
tion. However, for T <C T cn the superfluidity exponentially 



reduces the heat capacity. If the stellar crust was sufficiently 
colder than in Fig. ^, any neutron superfluid would reduce 
the contribution of free neutrons to negligibly small values, 
and the heat capacity would be determined by the ions and 
electrons. These effects are analogous in the core. The heat 
capacity in the superfluid core at T < T cn and T <C T cp 
would be determined by the electrons. 

The thermal conductivity in the core is taken as a sum 
of the conductivities of the electrons (Gnedin & Yakovlev 
1995) and neutrons (Baiko et al. 2000). The electron con- 
tribution usually dominates. The conductivity in the crust 
is assumed to be due to the electron scattering off atomic 
nuclei. In the outer crust the finite size of the proton charge 
distribution within a nucleus can be neglected. We use 
the recent results of Potekhin et al. (1999), which include 
the multiphonon processes in electron-nucleus scattering in 
Coulomb solid and incipient long-range nucleus-nucleus cor- 
relations in strongly coupled Coulomb liquid of atomic nu- 
clei. In the inner crust, we have performed original calcula- 
tions using the same formalism but taking into account the 
finite size of the proton core of atomic nuclei (finite-size nu- 
clei). Details of the calculations and the numerical fits are 
given in the Appendix. 

The results are illustrated in Fig. ^, which shows the 
density profiles of the electron thermal conductivity for 
T = 10 7 , 10 s and 10 9 K. The conductivity in the core is 
several orders of magnitude higher than in the crust since 
there are no such efficient electron scatterers in the core as 
atomic nuclei. In the crust we plot the thermal conductiv- 
ity for finite-size nuclei (solid lines) and for point-like nuclei 
(dashed lines). The finite-size effects are negligible near the 
neutron drip point but increase the conductivity at the crust 
base (p ~ 10 14 g cm -3 ) by a factor of 3 — 5. 

Another contribution to the thermal conductivity in 
the crust may come from the scattering of electrons off 
charged impurities, the randomly distributed nuclei of dif- 
ferent charge number. However, the most important temper- 
ature interval for neutron star cooling is ~ 10 8 — 10 9 K. We 
have verified that the thermal conductivity is almost unaf- 
fected by the impurities in this temperature interval (for a 
not very impure matter, Qi mp Si 1, where Qi mp is the impu- 
rity parameter defined in the Appendix) , although they can 
noticeably decrease the conductivity at lower T. Therefore, 
we neglect the effects of impurities in the present calcula- 
tions. 



3 THERMAL RELAXATION IN 

NON-SUPERFLUID NEUTRON STARS 

First, consider thermal relaxation in a young non-superfluid 
neutron star. The main features of the process are known 
from previous cooling simulations (Page & Applegate 1992; 
Lattimer et al. 1994; Page 1998a, 1998b; and references 
therein). 

Fig. |E| shows cooling models of neutron stars with dif- 
ferent masses. In the low-mass models, M < 1.44 Mq, the 
direct Urea process is forbidden. These stars follow the stan- 
dard cooling scenario and the cooling curves are almost in- 
dependent of M. The high-mass models go through the fast 
cooling scenario and demonstrate a spectacular drop of the 
surface temperature at the end of the thermal relaxation 
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Figure 6. Temperature profiles in the 1.3 Mq neutron star 
without supcrfluid effects. Numbers next to curves show stellar 
age. Contours arc at 0, 1CT 5 , 10~ 4 , 1(T 3 , 1CT 2 , 10" 1 , 1, 2, 5, 10, 
20, 50, 100, and 1000 yr. 
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Figure 8. Ratio of the integrated heat capacities in the crust 
and the core for the neutron star models of masses 1.1, 1.2, 
1.7 M Q . 
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Figure 9. Ratio of the neutrino luminosities in the crust and 
the core for the neutron star models of masses 1.1, 1.2, 1.7 
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epoch, t ~ 50 yr, due to the emergence of the cooling wave 
on the surface. The same, although much less pronounced, 
effect takes place in the case of slow cooling. 

Since the neutrino emissivity of the direct Urea process 
is several orders of magnitude larger than that of the mod- 
ified Urea processes, the fast cooling regime is established 
even if the central kernel, where the direct Urea process is 
allowed, occupies a small fraction of the stellar core (Page 
& Applegate 1992). For the high-mass stars (M > 1.44 M ) 
the cooling curves, again, depend weakly on the mass. The 
change of the slope of the cooling curves at t ~ 10 5 -10 6 
yr manifests the transition from the neutrino to the photon 
cooling stage. 

The surface temperature at the initial cooling stage (the 
first 100 years) is rather independent of the equation of state, 
stellar mass, or the core neutrino luminosity. The surface 
temperature is mainly determined by the physical properties 
of matter in the crust. The core and the crust are thermally 
decoupled, and the effective surface temperature does not 
reflect the thermal state of the stellar core. 

In contrast, the evolution of the central temperature, 
T(0, t), is drastically different for the slow and fast cooling 
scenarios. In all low-mass models, T(0,t) oc t^ 1 ^ 6 through- 
out the neutrino cooling era, t 10 5 yr, with a small offset in 
normalization. This follows from the simple power-law tem- 
perature dependence of the heat capacity [C v oc T) and the 
standard neutrino emissivity (Q v oc T 8 ). In the models with 
fast cooling, where the dominant neutrino process is the di- 
rect Urea (Qt, oc T 6 ), the scaling relation is T(0,t) oc t~ x ^ 
for the initial period t < 10 -2 yr. But then until t 5S 10 yr, 
the central temperature remains almost constant at 10 8 K 
as the heat flows from the warmer outer core, in which di- 



rect Urea process is prohibited, into the inner core. During 
the thermal relaxation epoch, 10 < t < 100 yr, the cen- 
tral temperature falls again by a factor of several. After the 
thermalization utill the end of the neutrino era, once again, 
T(0,t) oct- 1/4 . 

Figs. | and | illustrate the effects of thermal relaxation 
on the internal temperature profiles in the slow and fast 
cooling scenarios, respectively. Until the age of about 1 yr, 
the neutron star core, the inner and the outer crusts form 
almost independent thermal reservoirs. The region around 

4 x 10 11 g cm~ 3 , where free neutrons appear in the crust, 
seems to be the most effective at cooling owing to the pow- 
erful neutrino emission (see below). The outer crust cools 
to 10 9 K in less than a month, while the inner parts re- 
main much hotter. The core also cools independently but 
is unable to affect the inner crust due to the slow thermal 
conduction. During the first years the central kernel of the 
1.7 Mq model in Fig. ^ remains much colder than the outer 
core. This is because the kernel is cooled by the powerful 
direct Urea process and thermal conduction is still unable 
to establish thermal relaxation throughout the core. Almost 
full core relaxation is achieved in 10 years. 

After the first year, the crust temperature profiles of 
the slow and fast cooling scenarios start to differ. In the 
former, the temperature gradient between the core and the 
crust is slowly eroded, as the cooling wave from the center 
reaches the surface. In the latter, the temperature gradient 
continues to grow until it reaches a maximum at t ~ 10 yr. 
Then a huge amount of heat releases from the crust and 
leads to a spectacular drop of the surface temperature by 
an order of magnitude (which lowers the photon luminosity 
by four orders of magnitude). At t = 50 yr, the entire star 
is already isothermal. Note, that despite larger temperature 
gradients, thermal relaxation proceeds overall quicker in the 
fast cooling scenario. 

Prior to thermal relaxation, the contributions of the 
neutron star crust to the integrated heat capacity and neu- 
trino luminosity are significant (Figs. ^ and ^ . In the slow- 
cooling models, the heat capacity in the crust ranges from 
10% to 20% of that in the core, with the larger fraction in 
the low-mass models (where crusts occupy larger fraction 
of the volume). In the fast-cooling models, the ratio of the 
crust to core heat capacities reaches a maximum of 55% at 
t ~ 10 yr before dropping to under 10% after the relaxation. 
Similarly, the integrated neutrino luminosity of the crust is 
about 15%-40% of that of the core at t ~ 1 yr, and then it 
drops to a tiny fraction at later times. 

The importance of the individual neutrino mechanisms 
for the crust cooling varies at different epochs. First, for 
t < 10~ 2 yr in the fast cooling scenario or for t < 3 x 10 -3 
yr in the slow cooling scenario, the e~e + pair emission con- 
trols the crust temperature. As the temperature drops below 

5 x 10 9 K, this process quickly fades away. The next epoch is 
controlled by plasmon decay. It dominates for 10~ 2 < t 5o 10 
yr (fast cooling) or 3 x 10~ 3 < t 5s 10 yr (slow cool- 
ing). The last epoch of thermal relaxation lasts for the pe- 
riod 10 <; t < 100 yr (fast cooling) or 10 <, t < 1000 
yr (slow cooling), where either electron-nucleus or neutron- 
neutron bremsstrahlung is important. In fact, both neutrino 
processes give almost identical cooling curves in the ab- 
sence of superfluidity. However, free neutrons in the crust 
are thought to be in a superfluid state which strongly 
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suppresses nn bremsstrahlung. Therefore, electron-nucleus 
bremsstrahlung is likely to be the dominant neutrino mech- 
anism in this last epoch. 



4 RELAXATION TIME 

The duration of the thermal relaxation epoch is potentially 
interesting from the observational point of view. This prob- 
lem has been studied in a number of papers, with the most 
detailed and thorough work by Lattimer et al. (1994; also 
see references therein). Those authors considered thermal 
relaxation for the fast cooling and defined the relaxation 
time t w as the moment of the most negative slope of the 
cooling curve, lnT 3 (lnt), of a young neutron star. This is a 
typical time for the cooling wave to reach the surface. We 
will mainly focus on the case of rapid cooling, where the 
relaxation effects are more pronounced. 

According to Lattimer et al. (1994), the relaxation time 
of rapidly cooling neutron stars of various masses is deter- 
mined mainly by the crust thickness AJ? crus t and is given 
by a simple scaling relation 

t w »ah, a= (^^l) 2 (l-r 9 /i?)- 3 / 2 . (3) 

Here, ti is the normalized relaxation time which depends 
solely on the microscopic properties of matter, such as the 
thermal conductivity and heat capacity. In superfiuid neu- 
tron stars, t\ is sensitive to the magnitude and density de- 
pendence of the critical temperature of neutron superfluidity 
in the crust, as we discuss later. It is important that t\ ap- 
pears to be almost independent of the neutron star model, 
its mass M and radius R. We have verified that this scaling 
holds also for the slow-cooling non-superfluid models. 

For the non-superfluid stars with the core-crust inter- 
face placed at po/2, which is close to our value p cc = 1.5 X 
10 14 g cm" 3 , Lattimer et al. (1994) obtained h » 26 yr. Our 
rapidly cooling models show similar scaling, t\ — 28.4 ± 0.2 
yr, whereas in the slowly cooling models, t\ = 32.9 ± 1.2 yr. 

The dependence of t w on the thermal conductivity k 
and heat capacity C v follows from a simple estimate of the 
thermal relaxation time in a uniform slab of width 

t w ~ CJ 2 /k. (4) 

The proper width of a thin crust (A_R crus t <C R), tak- 
ing into a ccount th e effects of general relativity, is I — 
AR CTUSt /yJl — r g /R. This gives t w oc 1/(1 —r g /R) in equa- 
tion (||). An additional factor 1/ y/l — r g /R accounts for the 
gravitational dilation of time intervals. 

The parameter ti can be roughly estimated using equa- 
tion (0) with the values of the specific heat capacity and 
thermal conductivity in the crust near the crust-core inter- 
face, for instance, at p — 10 14 g cm" 3 . Fig. |^ shows that 
the relevant temperature at the interface is T ~ 2 x 10 8 K. 
At those temperature and density according to our physics 
input, C v = 1.9 x 10 19 erg cm -3 K _1 and k = 1.8 X 10 20 
erg cm" 1 s" 1 K" 1 . This gives ti ~ (1 km) 2 C„/K « 34 yr, 
in qualitative agreement with the model value. 

Figure [j^ shows the sensitivity of the relaxation time 
to test variations of physical properties of the neutron star 
crust in the fast cooling scenario of the 1.5 Mq star. The 
corresponding values of t w and ti are listed in Table B. They 



Table 2. Relaxation time t w and normalized time t\ for neutron 
stars with different crust models 



test crust model 


C (yr) 


*? (yr) 


t\ (yr) 


real model, no SF 


52.4 


28.8 


33.9 


no crust neutrinos 


253.5 


139.2 


134.9 


only plasmon decay neutrinos 


67.6 


37.1 


41.6 


only eZ neutrino bremsstr. 


58.3 


32.0 


34.5 


no neutron heat capacity 


15.3 


8.4 


6.7 


cond. for point-like nuclei 


131.8 


72.4 


102.3 


real model, weak crust SF 


20.2 


11.1 


3.3 


no Cooper neutrinos 


29.0 


15.9 


19.0 


weak core+crust SF 


22.3 


12.2 


25.7 


real model, strong crust SF 


15.0 


8.2 


6.7 


no Cooper neutrinos 


15.5 


8.5 


6.9 


strong core+crust SF 


10.7 


5.9 


5.8 


a for the 1.5 Mq model, with a 
b for the 1.3 Mq model, with a 


= 1.821 
= 2.875 
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Figure 10. Thermal relaxation for the 1.5 Mq model without su- 
perfiuid effects (see also Table ^). Solid line: the real cooling curve. 
Dotted lines: switched off neutrino emission from the crust (up- 
per) or infinite thermal conductivity at p > 10 10 g cm -3 (lower). 
The dashed curve C n = 0: removed neutron heat capacity in the 
crust. Another dashed curve: the thermal conductivity k in the 
crust is for point-like nuclei. Two other dashed lines: removed all 
neutrino mechanisms in the crust except either plasmon decay 
(pi) or electron-nucleus bremsstrahlung (eZ). 

may differ from the mean values (Table ^) within the error- 
bars. We present also the results for the superfiuid models 
discussed in the next section. 

Switching off the neutrino emission from the crust, while 
keeping the heat capacity, slows down the thermalization 
epoch by almost an order of magnitude, from 50 to 250 
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yr. Turning on the plasmon decay alone would lead to the 
cooling curve not very different from the real one, although 
the thermalization epoch would be delayed by about 30% 
(t w = 68 yr) . Turning off the plasmon decay and switching 
on the electron-nucleus bremsstrahlung instead would give a 
hotter neutron star before the relaxation but almost correct 
relaxation time, t w = 58 yr. Combined together, plasmon 
decay and electron-nucleus bremsstrahlung would reproduce 
accurately the real cooling curve. 

Restoring the neutrino emissivity in the crust but 
switching off the heat capacity of crustal neutrons (C„ = 0) 
we obtain much faster relaxation, which lasts about 15 years. 
This numerical experiment imitates the suppression of the 
neutron heat capacity by strong neutron superfluidity (dis- 
cussed in Sect. H). On the other hand, had we neglected the 
quantum suppression of the heat capacity of the nuclei be- 
low the Debye temperature, the latter heat capacity would 
have become important in older neutron stars, t <; 10 4 yr, 
strongly delaying the cooling. 

Finally, the relaxation time depends on the thermal con- 
ductivity in the inner crust. For instance, a neglect of the 
finite sizes of atomic nuclei in the electron-nucleus scatter- 
ing rate would lower the electron thermal conductivity at 
the crust base (p <; 10 13 g cm' 3 ) by a factor of 2-5 (cf Fig. 
^]). Using that, less realistic thermal conductivity we would 
have had much longer relaxation (about 130 years). If the 
thermal conductivity were infinite in the stellar interior (at 
p > pb), we would have obtained an isothermal cooling sce- 
nario. In this case a sharp drop of the surface temperature 
associated with the relaxation disappears. 

We have run additional cooling models with M = 
1.5 Mq, varying the heat capacity and thermal conductivity 
within the crust (at pb < p < p cc ) by a fixed factor 1/8, 1/4, 
1/2, 2, 4, 8. The results confirm that t w and ti are indeed 
quite sensitive to the variations of C v and k, in agreement 
with the qualitative estimate, equation (^), and the results 
of Lattimer et al. (1994). It is important that these vari- 
ations do not invalidate the scaling relation for the relax- 
ation time, equation (^[). Moreover, if C v is increased and n 
is decreased, the dependence of t\ on the values of the heat 
capacity and thermal conductivity is described by a simple 
scaling relation ti oc C v /k°' s . 

We have also done sensitivity tests of the relaxation 
time, analogous to those performed by Epstein et al. (1983). 
In each density region, 0.2 dex of log p, either the heat capac- 
ity, thermal conductivity, or neutrino emissivity have been 
changed by a factor 1/8, 1/4, 1/2, 2, 4, 8. Figure [H] shows 
the variations of t w when Q v and n are reduced by a factor 
of 8, and C v is enhanced by a factor of 8. The relaxation 
time depends mostly on the values of C v and k in the crust 
in the density range 10 13 — 1.5 x 10 14 g cm -3 near the crust- 
core interface, being rather insensitive to the variations of 
Q v . Variations of the physical parameters in the core af- 
fect the crustal relaxation much less strongly (at least for 
non-superfluid models) . The density range where t w is most 
affected by the variations of k is narrower than the den- 
sity range where it is affected by the variations of C„. The 
most important temperature range which influences t w is 
10 8 ^ T <, 10 9 K. The properties of matter in these "sen- 
sitivity strips" of p and T are very model-dependent. For 
instance, the nuclei may be strongly non-spherical (rods, 
plates, etc.) at p > 10 14 g cm -3 , which is not included in 
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Figure 11. Sensitivity of the crust relaxation time to the vari- 
ations of the heat capacity, thermal conductivity, and neutrino 
emissivity in various density regions of the 1.5 Mq neutron star 
model with the non-superfluid core. In each 0.2 dex of logp, Q v 
and k are reduced by a factor of 8 (solid lines, no superfluidity), 
and C v is enhanced by a factor of 8 without superfluid effects 
(dashes) and with weak superfluidity (dots). 



our calculations. Note that the thermal conductivity has not 
been calculated so far for the phase of non-spherical nuclei. 

However, if G v were noticeably lower or k noticeably 
higher than in our basic non-superfluid models, the decrease 
of the crust relaxation time would saturate at t m ~ 13 yr. 
This is the time it takes the inner core with the direct Urea 
emission to equilibrate thermally with the outer core (cf Fig. 
Q). More generally, this is a core relaxation time t COIC , which 
is almost independent of the parameters of the crust. It can 
be estimated using the same formalism of heat diffusion, 
equation (Q) , through a slab of material between the direct- 
Urca-allowed kernel and the boundary of the core, with 
I = Rcora — Rd- The relativistic factors appear as (e~* cc ) 3 , 
where $ cc is the metric function at the core-crust interface. 

We have calculated the core relaxation time by setting 
the crust conductivity very high for the fast cooling models 
of non-superfluid neutron stars with M — 1.5, 1.6, 1.7 Mq. 
The least-squares fit gives 



f Rcora — Rd \ 

\ 10 km J 



(5) 



with — (9.1 ± 0.8) yr. Thus, the fast cooling models may 
have two distinct relaxation times, in the core and in the 
crust, and the latter is typically longer, at least for non- 
superfluid models. 

The normalized crustal relaxation time ti shows very 
small variations with the neutron star mass, ±0.2 yr for fast 
cooling and ±1.2 yr for slow cooling (see Table ^|). Even 
combining the models with all masses we find that the scal- 
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Figure 12. Temperature profiles in the interior of the 1.3 Mq 
model with (solid lines) and without (dots) weak crust superflu- 
idity of free neutrons. The core is not superfluid. Numbers next 
to curves show the stellar age. The contours are at 0, 10 — 4 , 10 -2 , 
1, 100, and 1000 yr. 
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Figure 13. Temperature profiles (solid lines) in the interior of 
the 1.5 Mq model with strong superfluidity both in the crust 
and the core. Numbers next to curves show the stellar age. The 
contours are at 0, 10 — 4 , 10 -2 , 1, and 20 yr. Dotted lines show 
the temperature profiles of the non-superfluid star. 



ing relation, eq. (|3|), holds remarkably well, ii(all) = 31 ± 2 
yr. Therefore, the cooling of the crust prior to thermal re- 
laxation is indeed insensitive to the thermal evolution of 
the core. The slow cooling makes the relaxation much less 
pronounced but it does not change its basic features. This 
property owes to the very high thermal conductivity in the 
core (Fig. ^) , which makes the core relaxation time typically 
much shorter than the crust one. Accordingly, the cooling 
curves do not depend on the exact values of the conductivity 
in the core as long as it is much higher than in the crust. 



5 EFFECTS OF SUPERFLUIDITY 

Free neutrons in the crust and both neutrons and protons in 
the core of a neutron star are likely to be in the superfluid 
state. We im plement the effects of superfluidity as discussed 
in Sect. 



2.2 . It turns out that superfluidity in the crust affects 



the cooling curves at the thermal relaxation stage, while su- 
perfluidity in the core affects cooling at later stages. First, 
consider the effects of neutron superfluidity in the crust. 
If the temperature T falls much below the critical tempera- 
ture Ten, the superfluidity reduces strongly the neutron heat 
capacity and nn neutrino bremsstrahlung. While the latter 
is compensated by electron-nucleus bremsstrahlung, the for- 
mer effect leads to a faster thermal relaxation. In addition, a 
new neutrino process is allowed in the superfluid state, the 
neutrino emission due to Cooper pairing of free neutrons. 
It further accelerates cooling and thermal relaxation of the 
crust. 

Figure ^ illustrates how the weak neutron superfluidity 
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Figure 14. Superfluid effects in the crust of the 1.5 Mq neutron 
star model with the non-superfluid core. Dashed line is for the 
case of weak neutron superfluidity, while dash-and-dot line is for 
the case of strong superfluidity. Dotted line is obtained neglecting 
the Cooper-pair neutrino emission for the weak-supcrfluid model. 
Solid line is the cooling curve for the non-superfluid crust. 
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in the crust carves out the temperature profiles in the stan- 
dard cooling scenario, M = 1.3 Mq. For the first 10 -4 yr, 
when the temperature is above T cn , the thermal structure 
is identical to that of the non-superfluid model. Later, the 
region where the critical temperatures are the highest cools 
much faster than the neighboring layers. The acceleration of 
cooling is mainly due to the reduction of the heat capacity 
and switching on the Cooper-pairing neutrino emission. A 
sequence of points in Fig. ^ in which solid lines start to 
deviate from the dotted ones reproduces the density profile 
of T cn shown in Fig. [lj As the temperature falls further, 
wider density regions become affected, producing shells of 
cool matter surrounded by hotter layers on both sides. Af- 
ter the cooling wave from the core reaches the outer crust, 
the star settles into almost the same isothermal state as the 
non-superfluid model, but faster. 

The effects are much stronger if the superfluidity is al- 
lowed for in the stellar core. Figure [l^ shows the combined 
effect of the strong superfluidity in the core and crust. Both 
neutrons and protons are superfluid in the core of this 1.5 
Mq neutron star. In addition to the trough in the crust lay- 
ers, the core develops a complex structure. All the sources 
of neutrino emission and the nucleon heat capacity in the 
core are affected by the superfluidity, while the electron heat 
capacity is not and becomes dominant in highly superfluid 
cores. As soon as the temperature drops significantly below 
the critical temperatures T cn or T cp , Cooper-pairing neu- 
trino emission comes into play. It starts in the inner part of 
the core and drives the temperature down. At t ~ 10~ 2 yr 
that region is even cooler than in the non-superfluid model, 
while the other parts of the core are slightly hotter. By the 
age of 1 yr, this cool region includes all of the core except 
the inner kernel. As a result, thermal relaxation proceeds on 
a shorter timescale, and at t 20 yr the stellar interior is 
isothermal. Thus, for the large assumed values of T cn and 
T cp , the neutrino emission due to Cooper pairing becomes 
so strong that, instead of slowing down, the presence of the 
core superfluidity accelerates the cooling. 

Figure [w| demonstrates the aforementioned effects on 
the thermal relaxation time in the 1.5 Mq neutron star. 
The stellar core is assumed to be non-superfluid. The ther- 
mal relaxation stage occurs about 2.6 times faster in the 
case of weak crust superfluidity and about 3.5 times faster 
in the case of strong crust superfluidity, compared to the 
non-superfluid crust. And also, while the inclusion of the 
neutrino emission due to Cooper pairing leads to faster cool- 
ing for t < 20 yr, most of the accelerating effect is due to 
the reduction of the crust heat capacity. Therefore, the re- 
laxation time is greatly reduced by the effect of superfluidity 
on heat capacity of free neutrons in the stellar crust. 

Still, the relaxation time of the superfluid crust satis- 
fies the same scaling relation, equation (^), as for the non- 
superfluid crust. Table |^ summarizes the results, taken as 
the mean of M — 1.5, 1.6, 1.7 Mq models for fast cooling, 
and M = 1.1, 1.2, 1.3, 1.4 Mq models for slow cooling. In 
case of the fast-cooling models, crust superfluidity reduces 
the relaxation time by up to a factor of four; t\ — 7.0 ± 0.9 
yr for the strong superfluidity. In case of the slow-cooling 
models with strong crust superfluidity the effect is similar, 
but with the weak superfluidity it is much different, almost 
a factor of 10. The inclusion of superfluidity in the core does 
not affect the relaxation time significantly. The values of t\ 



Table 3. Normalized relaxation times in the crust t\, eq. (J3J) 
in the core t2, eq. (H) 



and 





fast cooling 


slow cooling 


crust ti (yr) 






no SF 


28.4 ±0.2 


32.9 ± 1.2 


weak crust SF 


10.3 ±0.6 


3.4 ±0.3 


weak core+crust SF 


11.5 ±0.5 


25.0 ±0.4 


strong crust SF 


7.0 ±0.9 


6.8 ± 1.0 


strong core+crust SF 


6.2 ±0.4 


5.7 ±0.2 


core ti (yr) 






no SF 


9.1 ±0.8 




weak core SF 


11.2 ±0.5 




strong core SF 


3.6 ±0.2 





change by about ~ 1 yr and are within the error-bars of the 
results for the non-superfluid cores. The notable exception 
again is the weak superfluidity in the slow-cooling models. 

Superfluidity in the core changes the core relaxation 
time, although it still scales as equation Weak core su- 
perfluidity only increases the core heat capacity and delays 
the relaxation by a few years, ti = 11.2 ± 0.5 yr, while the 
strong superfluidity suppresses the heat capacity and gives 
much faster relaxation, by a factor 2.5. 

The large difference of the values of t\ for the slowly and 
rapidly cooling stars with weakly-superfluid crusts indicates 
that, generally, superfluidity may violate the simple scaling 
relations for the relaxation time as well as the strict thermal 
decoupling of the crust and the core prior to the relaxation. 
This happens because superfluidity makes the heat capacity 
(and other properties of matter) a strong function of density. 

The shortening of the thermal relaxation phase in a 
rapidly cooling star, due to the superfluid reduction of the 
crustal heat capacity, was emphasized by Lattimer et al. 
(1994). Using a model of superfluidity, which is closer to our 
model of strong superfluidity in the crust, Lattimer et al. 
find that the relaxation time becomes three times as short, 
t w — 8.4 ± 2.0 yr, in qualitative agreement with our results. 
Our calculations indicate that the effect is sensitive to the 
model of neutron superfluidity in the crust and thus, it can 
be used to test such models. 

When thermal relaxation is over, the cooling is mainly 
regulated by the neutrino luminosity and heat capacity of 
the core. The processes in the crust cease to play a significant 
role except for the very low-mass neutron stars with large 
crusts. We find, however, that for some superfluid models 
the neutrino luminosity of the crust may affect the cooling 
for a short period of time during the transition from the 
neutrino cooling era to the photon era. For instance, cooling 
of superfluid neutron stars after the thermal relaxation is 
discussed in detail by Yakovlev et al. (1999). 



6 CONCLUSIONS 

We have studied thermal evolution of young neutron stars 
using a new cooling code and updated physics input. The 
effective surface temperature of isolated neutron stars in the 
first 100 yr is determined by the properties of the crust. Af- 
ter the cooling wave from the core reaches the surface, the 
effective temperature drops from 250 eV to 50 eV or lower. 
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We confirm the conclusion of Lattimer et al. (1994) that the 
duration of the relaxation epoch in a fast cooling scenario 
depends sensitively on the heat capacity and thermal con- 
ductivity of the crust. The relaxation time scales with the 
size of the crust, eq. (|^), and the normalization is t\ = 28 
yr without superfluid effects or t\ — 7 yr with strong su- 
pcrfluid effects. We find that the same conclusion holds in 
the case of slow cooling, in which thermal relaxation is less 
pronounced. 

We also investigate the effects of various neutrino emis- 
sion mechanisms in the crust, superfluidity of free neutrons, 
heat capacity, and thermal conductivity on the relaxation 
time (as summarized in Table ^| and Fig. [y]). The relax- 
ation time is most sensitive to the variations of k and C v in 
the density range 10 13 -1.5xl0 14 g cm 3 near the crust-core 
interface. 

Young cooling neutron stars can serve as astrophysical 
laboratories of matter at subnuclear densities. The sensitive 
dependence of the relaxation time on the microscopic prop- 
erties of matter in the deep inner crust provides a possibility 
to study these properties by observing the emergence of the 
cooling wave. In order to realize this method, one needs to 
detect thermal emission from a very young neutron star in 
the range from 50 to 250 eV. Such stars have not been found 
so far, but may be detected in young supernova remnants 
by Chandra and XMM. 

To improve the cooling theory of young neutron stars, 
it would be desirable to continue theoretical studies of mi- 
croscopic properties of matter in the inner crust, in par- 
ticular, the possible unusual phases of nonspherical nuclei 
at p > 10 14 g cm -3 , as well as the thermal conductivity, 
heat capacity, and neutrino emission of such matter. It is 
also important to refine our knowledge of matter contain- 
ing spherical nuclei at p £ 10 14 g cm -3 . In this paper we 
present (Appendix) and use the results of new calculations 
of the electron thermal conductivity due to the scattering 
of electrons off spherical atomic nuclei, taking into account 
the effects of finite sizes of proton nucleus cores. 
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APPENDIX A: ELECTRON CONDUCTIVITY 
OF THE INNER CRUST 

In the inner crust, the electrons are strongly degenerate. The 
heat and charge are transported mainly by the electrons with 
energies in a narrow thermal band near the Fermi level ep. 
The coefficients of the electrical and thermal conductivities 
can be written as 



(Al) 



rn* 3me 

where m* = £f/c 2 , n c is the electron number density, fee is 
the Boltzmann constant, and r CT and r K are the effective re- 
laxation times. In the presence of several electron scattering 



mechanisms, the cumulative relaxation time is determined 
by the Matthiessen rule (e.g., Ziman 1960). In the envelopes 
of neutron stars, one traditionally considers the electron- 
ion (ei), electron-electron, and electron-impurity scatterings. 
We focus on the main mechanism, the ei scattering. 

A modern theoretical treatment of the conductivi- 
ties due to the ei scattering in strongly coupled Coulomb 
plasmas of ions (atomic nuclei) has been proposed by 
Baiko et al. (1998) and applied to the outer envelopes 
of neutron stars by Potekhin et al. (1999). In the lat- 
ter work, practical formulae for a and k incorporating 
all three scattering mechanisms have been derived; the 
corresponding Fortran code is availab le electronically at 
http: / / www.ioffe.rssi.ru / astro / conduct / . 

For the inner crust these expressions must be modified. 
First, the atomic nuclei cannot be considered as point-like. 
Second, the ei scattering in the inner crust at temperatures 
much below 10 8 K changes its character: so called umklapp 
processes cease to dominate and the normal processes (with 
electron momentum transfer within one Brillouin zone) be- 
come more important. 



Al Effects of finite size of ions 



The effective ei relaxation time is 
3ivh 



4Ze F Q f 2 A ei (e P )' 



(A2) 



where at = e 2 /(he) is the fine structure constant, and A e i 
is the Coulomb logarithm. For a classical plasma of ions in 
the Born approximation (e.g., Potekhin et al. 1999), 



A, 



dq q 3 \(f>(q)\ 2 S(q) 



1 2k F 



(A3) 



4>(q) 



Hi) 

q 2 e(q) : 



where &f = (3n 2 n e ) 1 '' 3 is the electron Fermi wave number, 
/3 r = vf/c, vf = hkF/mt is the Fermi velocity, e(q) is the 
static dielectric function, S(q) is the static structure factor 
of the ions (more precisely, its inelastic part - see Baiko et 
al. 1998), and F(q) is the form factor of the ions. The latter 
three functions allow for the electron polarization, ion-ion 
correlations, and the finite size of the ions, respectively. 

The quantization of ionic motion b ecomes important at 
T < T p , where T p = h ^^TTmZ 2 e 2 /M /k B is the ion plasma 
temperature, rii = n a /Z is the ion number density, M — 
A m u is the ion mass, and m u is the atomic mass unit. In this 



case, in equation (A3) one should use an effective structure 



factor, different for a and k. These effective factors have been 
derived by Baiko et al. (1998) and fitted by Potekhin et al. 
(1999) for the case where the umklapp processes dominate. 



A 1.1 Uniform charge distribution in atomic nuclei 

At p <C p cc a good approximation for the charge density 
within a nucleus is a step function. The corresponding form 
factor is 

3 

F(q) = ^- [sin(qT nuc ) - gr nuc cos(gr nU c)] , (A4) 

\qr nuc J' 

where r nuc is the proton core radius. 
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In general, A e i depends on p, T, Z, A, r nuc , and it ap- 
pears to be almost independent of the type of Coulomb 
lattice. It is convenient, however, to introduce the dimen- 
sionless parameters: the electron relativity parameter x T = 
hkp/(m e c); the ion size parameter a; lluc = r nuc /a,i, where 
ai = (47rni/3)~ 1//3 is the ion-sphere radius; the Coulomb 
coupling parameter of ions T = Z 2 e 2 /(aiksT); the inverse 
ion quantum parameter t p = T/T p ; and the electron screen- 
ing parameter s c = k^ F /(2k F ) 2 = ctf/(-Kf3 r ), where fcxF is 
the electron screening wave number. In accordance with the 
previous results, it is also convenient to introduce an auxil- 
iary parameter sd = (2k F ?-d)~ 2 , where td = ai/s/ST is the 
Debye screening length for the ideal plasma of ions. Finally, 
the basic parameters that characterize the Coulomb lattice 
are the normalized first- and second-negative moments of 
the phonon frequency, it_i ~ 3 and it_2 ~ 13. 

We have calculated numerically the effective Coulomb 
logarithms for about a hundred pairs (p,T), varying 

over orders of magnitude in the domain of strongly degener- 
ate electrons and strongly coupled ions, for eight ion species 
from Z = 12 to Z — 60, and for five values of x nuc from to 
0.4. We have included the non-Born corrections in the same 
manner as Yakovlev (1987). The results can be fitted using 
the effective screening function 



-(1 + sw) e aw [Ei (sw) - Ei (sw + w) 



e3 (l)\ 2 = 



(q 2 +qi? 



Go.k D 



(A5) 



instead of \<f>(q)\ 2 S(q) in equation (Ai). Here, the first term 
corresponds to the Debye screening with the effective screen- 
ing wave number q 3 , the factor in square brackets corrects 
for the ion correlations, and the functions G and D describe 
ionic quantum effects, as in Potekhin et al. (1999). An addi- 

_ 2 

tional factor e siq plays a role of the effective form factor. 
The numerical values of A obtained from the accurate theory 
are reproduced (within several percent in the most impor- 
tant p—T range) if we use the effective screening function 
(|A5j) with the following parameters: 



q s /(2k F ) = (si + s c 
Pz = ira{Z(3 T , 



-Pz 



Si = s D (1 + 0.06 T) e 



-vT. 



w = (2fe P ) 2 so = (u-s/sd) (1 + 0z/3); 



Wl 



so 



(2k F 



14.73 xt 



G a = (1 + 0.0361 z~ lf yt p ) 
0.0105 t p 



G k — Go ~\~ 



(t 2 + 0.0081) 3 / 2 



Z/13) (l+p z /3), 
1/2 (1 + 0.122/31), 

[l+tf(3 z ] (l-Z- 



D 



x (1 + a; nuc v2Z), 
exp -0.42 yJx r /{AZ)u--i exp(-9.1t p ) 



(A6) 
(A7) 

(A8) 
(A9) 

(A10) 
(AH) 



If inuc = 0, these equations reduce to those presented in 
Potekhin et al. (1999). 

With the effective screening function in the form 
we can integrate equation (A3) analytically and obtain 



= [Ao(s,w + tw)-Ao(s,™i)] G a , K D, (A12) 
where 



Ao(s,u>) = Ai(s,to) — /3 t A2(s, w), 



(A13) 



2Ai(s,w) = In ■ 



+ 



s + 1 



2 A 2 (s, w 



1+w 



(1 



2s In 



(A14) 
s + 1 



w s + 1 ' s 

+s (2 + sw) e sw [Ei (sw) - Ei (sw + w)] , (A15) 

and Ei (a;) = y~ 1 e~ v dy is the standard e xpon ential in- 
tegral. Note that using equations (A14) and QAlq) directly 
may result in large numerical round-off errors in the limit- 
ing cases s <C 1, w <C 1, or w 1. In these cases, explicit 
asymptotic expressions of Potekhin et al. (1999) should be 
used. 



A 1.2 Realistic charge distribution 

Near the bottom of the crust, the boundary of the proton 
core in nuclei becomes fuzzy and the above results should be 
modified. Oyamatsu (1993) approximated the proton charge 
distribution by a function proportional to (1 — r/r m ) b (at 
r < r m ), where the power index b controls the "sharpness" 
of the charge profile. The parameters b and r m have been 
described by the simple functions of mass density in the 
smooth composition model of crust matter (Kaminker et 
al. 1999). Using that model we have calculated the electri- 
cal and thermal conductivities for p = 10 9 — 10 14 ' 1 gcm~ 3 
and T = 10 7 — 10 9 K, and compared the results with the 



fitting formulae (A12)— ( A15). We have found that the elec- 
trical conductivity is reproduced (within w 5-10%), if we 
define x nuc in equation (A.8) so as to reproduce the same 



mean-square proton-core radius as in the approximation of 
uniform charge distribution: 



1 - 15/(5 + b) + 15/(5 + 2b) - 5/(5 + 36) 
1-9/(3 + b) +9/(3 + 26) -1/(1 +6) 



(A16) 



In equation (A1C) which determines the thermal conductiv- 
ity at T < T p , one should additionally multiply a; nuc by 
a factor 6/(0.5 + 6). The maximum fit error of k does not 
exceed 13%. 



A2 Normal processes at very low temperatures 

Near the boundaries of the Brillouin zones the dispersion re- 
lation of electrons differs from the free-electron case, and at 
the boundaries the electron energy spectrum contains gaps. 
The gaps Ae can be estimated in the weak coupling approx- 
imation (e.g., Kittel 1986) as Ae ~ <j>(k F ) — AnxZe 2 n\k F 2 = 
4e 2 /(37rfc P ). The effect of gaps is most significant if the de- 
viation of the electron momentum from the intersection line 
between the Fermi surface and the Brillouin zone boundary 
does not exceed Ak ~ Ae/(?iwp) ~ ^ (ai/(3 r ) k F <§C fcp. 
However, with decreasing temperature the strips of the 
Fermi surface, between which the umklapps proceed effec- 
tively, become narrower and closer to these intersection lines. 
When the widths of the strips, ~ t p (67r 2 ni) 1//3 , become 
smaller than Ak, the umklapp processes are frozen out and 
the normal processes prevail. The above estimates indicate 
that this happens when the temperature falls below 



T u ~ T p Z 1/3 a c /(3&). 



(A17) 



In this case, the formalism used in Sec. Al becomes 



invalid. The asymptotes for the ei scattering rates at T < T u 
have been obtained by Raikh & Yakovlev (1982). In our 
notation, they yield: 
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Figure Al. Density dependence of the electron thermal con- 
ductivity at T = 10 7 , 10 s and 10 9 K in the neutron star crust. 
Solid and dot-dashed lines: ground state composition (Negele & 
Vautherin 1973), rectangular profile of nuclear charge (solid lines: 
pure crystal or liquid; dot-dashed line: including impurities; bot- 
tom of the hatched region: neglecting freezing of umklapp pro- 
cesses at T = 10 7 K); dashed lines: smooth-composition model, 
realistic profile of nuclear charge, no impurities. 



K 
A!: 



a,Q x 



1/2 



where 



180 



/ 3 m e \ 
\irai m u ) 



(4/3) (a f //3 r )*p 
t 3 



(A18) 



1/2 



C(5) « 50, 



and £(5) = 1.0369 is the value of the Riemann zeta function. 

Now we interpolate between the high -tem perature 
Coulomb logarithm A^^ igh given by equation (A12) and the 
low-temperature asymptote A^ 1 ? given by equation (A18): 



Al 



= AX Ugh exp(-T u /T) 

+K&m [1 - exp(-T u /T)] . 



(A19) 



Figure Al shows the thermal conductivities for three 



values of the temperature. In the neutron star envelope, the 
plasma is in the solid phase except for the segment to the 
left of the asterisk at the highest T = 10 9 K. At the low- 
est T — 10 7 K, the change in the scattering from umklapp 
to normal processes leads to a significant increase in the 
conductivity shown by the hatched region. However, this 
increase can be compensated by the scattering off impuri- 
ties. For example, the dashed line shows the conductivity 
for Qi m p = {(2 imp - Z) 2 n imp ) /m = 0.2, where n imp and 
Zimp are the impurity number density and charge number, 
respectively. For comparison, dashed lines show the thermal 
conductivity in the smooth-composition model with smooth 
proton charge profile. 
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